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Abstract 



It is shown that the absolute value of the persistent current in a system 

with toroidal geometry is rigorously less than or equal to ehN/Airmr 2 , where 

N is the number of electrons, and r$ 2 = (r^ 2 ) is the equilibrium average of 

the inverse of the square of the distance of an electron from an axis threading 

the torus. This result is valid in three and two dimensions for arbitrary 

interactions, impurity potentials, and magnetic fields. 
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The phenomenon of persistent currents occurs when an electronic system is placed in 
a magnetic field: at thermal equilibrium, an electric current flows without dissipation of 
energy. This effect is usually studied in systems which are topologically equivalent to a 
torus, for example a metal ring, or a hollow cylinder. The interesting quantity is the flow of 
current through a cross section of the torus. 

The persistent current exhibits a variety of behaviors, depending on both the magnetic 
field and the geometric parameters of the system. A first example is that of a thin metal 
ring, i.e., a ring whose thickness is much smaller than its radius. Mesoscopic versions of this 
system have received particular attention in the past few years. The experiments [|IJ] were 
done in the Aharonov-Bohm configuration, in which a weak magnetic flux threads the ring, 
without significantly affecting the electrons orbits. In this case, the main physical effect 
is the "twisting" of the boundary conditions on the electron wave function, leading to a 
persistent current which is a periodic function of the threading flux, with period $o — hc/e 
0. The typical magnitude of the current is ev/L, where v is a characteristic velocity of 
propagation of an electron, and L is the length of the ring. The current vanishes if the ring 
is made larger and larger (L — > oo), so this effect is a purely mesoscopic one. 

A different behavior is obtained in a two-dimensional ring, with inner and outer radii 
Ri and Ri [R\ < R2), in the presence a strong perpendicular magnetic field B such that 
the magnetic length I = (hc/eB) 1 ^ 2 « R 2 — R\. This model has been studied by several 
authors The magnetic field induces currents flowing in opposite directions at the inner 
and the outer edges of the ring. If Ri and R2 are macroscopically different, the edge currents 
do not cancel each other exactly, and one is left with a net current that fluctuates violently 
as a function of of electron number when the Fermi level is in a gap between two Landau 
levels |3J]. A typical value of the order of a fraction of eu c (u c = eB/mc is the cyclotron 
frequency) has been reported in a numerical study Q. 

In the special limit R\ — > 0, the ring becomes a "punctured disk", and then it has been 
found || (neglecting disorder and interactions) that the net current is quantized in integral 
multiples of euj c /Air when the chemical potential is pinned to one of the Landau levels in the 
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bulk. 



In view of the diversity exemplified above, it is remarkable what we show in this paper, 
that there exists a rigorous upper bound to the persistent current of a system of arbitrary 
size and shape, provided that it is topologically equivalent to a torus. We do not need to 
assume any symmetry, and we do not put any constraints on the nature of the magnetic 
field, not even that it be uniform. In essence, our derivation of the upper bound, is a 
sharpening of the argument presented by Bohm ||, following a suggestion by Bloch, to 
prove that a macroscopic one- dimensional ring cannot carry a finite circulating current, at 
thermal equilibrium. Assuming that the upper bound is a good estimate of the maximum 
value of the persistent current that can be reached in a given system with an appropriate 
magnetic field, this result enables us to easily understand the large difference in order of 
magnitude and geometric dependence of persistent currents in, for example, thin rings and 
punctured disks. Also, the rigorous upper bound can be useful as a test of the validity of 
approximate theories of persistent currents. 

Let us consider a system of electrons confined within a body of toroidal topology, such 
as the one shown in Fig. 1. No symmetry is assumed. Let us choose an axis - the "z" axis 
- which threads the body, but is otherwise arbitrary. The question of the optimal choice 
of the z axis will be addressed later. Each half-plane emerging from the z axis at an angle 
(j> (0 < 4> < 27r) cuts a two-dimensional cross section in the body (see Fig. 1) @. The 
persistent current is defined as the flux of the current density j(f) through S^. By virtue of 
the continuity equation V • j(r) = 0, this flux is independent of the choice of (f). 

We introduce the standard cylindrical coordinates = (rj,Zj,0j) to characterize 
the position of the i-th electron in the body. The corresponding momenta are pi = 
—ih(d/dri,d/dzi,r^ d/d(f)i). The Hamiltonian of the system is 



where V(r) and A(r) are arbitrary scalar and vector potentials. The exact eigenfunctions 




(1) 




r N , zn, ) of H, with eigenvalue E n , are completely antisymmetric, and 
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vanish very rapidly (exponentially) outside the boundaries of the body. The average persis- 
tent current, at thermodynamic equilibrium, at temperature T = l/k B j3 is given by 

I = ^E^^nl / j*(?)drdz\1> n ), (2) 

Zj n J Sty) 

where Z is the partition function, and 

^(r) = ~ £{n^(r - fj) + 5(r - f,)n^}, (3) 

is the azimuthal component of the current density operator. The azimuthal component of 
the kinetic momentum operator is defined as 

- ih d e . ,_ s , d . 

As we have already remarked, the integral, in Eq. ([|) is independent of the angle character- 
izing the cross section S(<p) in the (r, z) plane. Using this fact, together with the definition 
of the current density operator, Eq. ([|), it is easy to verify that 

(5) 

where 

i = — ^v^, (6) 

% 1 

is the current operator, and (...) denotes the usual thermal equilibrium average, with Hamil- 
tonian H. 

The upper bound to / is derived as follows. Consider a gauge tranformation 

Ui = exp [il^fa], (7) 

i 

where I must be an integer in order to preserve the single-valuedness of the wave functions. 
The Hamiltonian is transformed to 

- — - Ih* h 2 l 2 1 

H' = U^Uf 1 = H--I + —T,± (8) 



Let F and F' denote the free energies associated with the Hamiltonians H and H' respec- 
tively. They satisfy the well known inequality M 



F' <F + (H'-H). 



(9) 



But, in this case, F' = F because H and H' are related by a unitary transformation. 
Therefore, combining Eqs. (pf) and (R), we obtain the inequality 



h z l 



2/2 
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Ih 



which must be satisfied for an arbitrary value of the integer I. Let us define 

- = -(Y-) 



(10) 



(11) 



and 



ehN 



(12) 



Then Eq. (flOf) takes the form 



I 2 -I— > 0. 



(13) 



It is now easy to verify that, for |/| > I m , Eq. ( |T3"D is always violated either by I — 1 or by 
I = —1. On the other hand, for |/| < I m , Eq. (^3[ ) is satisfied by all integral values of I. We 
conclude that 



l/l < /„ 



(14) 



the result of this paper. 

The above argument is also straightforwardly applicable to a strictly two dimensional 
ring geometry. The arbitrary axis threading the toroidal body is replaced by an arbitrary 
point within the hole, and the cylindrical coordinates relative to that axis are replaced by 
planar polar coordinates relative to this point. The cross section becomes a segment (or a 
collection of segments) of straight line, and the z coordinate is suppressed. The final result 



is still given by Eq. (0). Incidentally, the fact that the threading axis does not "touch" the 
toroidal body guarantees that I/tq is a finite number, because the wavefunction vanishes 
exponentially in the regions where r% — > 0. 

Equation ( |T4| ) is still dependent on the arbitrary choice of the reference axis. We can 
completely remove this arbitrariness, by choosing, for any given system, the axis that yields 
the most stringent inequality, i.e. the smallest value of 1/tq. Thus, our original definition of 
r in Eq. (|TT|) is replaced by 

1 11 

= min { )}, (15) 

where the minimum is calculated with respect to the set of all possible axes threading the 
toroidal body. Qualitatively, this is the axis which, on the average, remains as far as possible 
from the body. Notice that 1/tq < l/r^ in , where r min is the minimum distance from the 
optimal axis to the body. Since this latter quantity is independent of the magnetic field, 
Eq. ( p^D implies that the current as a function of magnetic field cannot increase monotoni- 
cally, but must either oscillate or tend to a constant (at fixed particle number). 

It is worth noting that the upper bound on the persistent current can be expressed in 
terms of a moment the electronic density distribution: 

eh f n(x) 

where r is the distance of point x from the coordinate axis. Thus, the calculation of I m in 
specific cases requires only a knowledge of the electronic density distribution - a relatively 
"gross" property of a many-electron system. Furthermore, from the form of eq. ([16]), it is 
clear that the value of I m is not sensitive to fine details of the density distribution, especially 
in bodies of macroscopic to mesoscopic size. Essentially, I m is a geometric parameter. 
In three dimensional macroscopic bodies, we expect that approximating the density as a 
constant within the geometric boundary of the body and zero outside, should lead to an 
excellent approximation for I m . A better estimate of the density distribution can be obtained 
by solving the Thomas-Fermi or the Kohn-Sham equations. 
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Let us now consider a few concrete applications of Eq. (p~4|) . Suppose that the body 
is a circular ring of radius R and negligible thickness ("thin ring"). Then the optimal axis 
coincides with the axis of the ring, and r$ ~ R. The upper bound to the current is in this 
case 

(17) 

where L = 2ttR is the circumference of the ring. This tends to for L — > oo, so there is 
no current flowing in a macroscopic thin ring. It must be noted that Eq. ( |17|) is rigorous 
for a geometrically thin ring, i.e., it incorporates the effect of electron-electron and impurity 
scattering exactly. The ring may still be three-dimensional as far as the electronic wave 
function is concerned. No assumption on the form of the density distribution has been used. 
In the special case of a non-disordered one-channel ring, containing non-interacting electrons 
with a one- dimensional Fermi velocity vp, the upper bound reduces to 2cvf/L (the factor 
2 arising from spin degeneracy). This value is attained exactly at T = 0, when the ring is 
threaded by a Aharonov-Bohm flux $ = hc/2e (odd N) or $ = + (even N). 

As a second example, consider a two dimensional disk (outer radius R2) with a central 
hole of radius Ri, in a strong magnetic field, such that / << R 2 — R\- Once again, the 
optimal axis coincides with the axis of the disk. Since the electronic areal density hq = N/A 
(A is the area of the ring between the inner and the outer edge) is essentially uniform in 
the bulk of the system, we can calculate 1/Vq by a taking a uniform average of 1/r 2 , which 
yields the result 1/rjj = 27rln (R 2 /Ri)/A. Substituting this in Eq. ([12"D , we obtain 

eh , , Ro . 

The behavior of the upper bound depends only on the aspect ratio R 2 /Ri. If R\ ~ R2 to 
an accuracy much smaller than Ri or R 2 , we simply recover the thin ring geometry, and the 
persistent current vanishes. In the limit R\ — *■ (or R 2 — > 00), we obtain the "punctured 
disk" geometry. The upper bound diverges logarithmically, allowing a large persistent cur- 
rent. (The logarithmic divergence is of course cut off when Ri becomes comparable to the 



size of the edge region, in which the approximation of constant density is no longer valid.) 
In all other cases, the logarithmic factor is a constant of order 1, and the upper bound to 
the current can be written as I m ~ eh/2irma 2 , where a is the average distance between the 
electrons. 

In the interesting case that the filling factor v = 2irl 2 no, rather than the areal density, is 
kept constant in the bulk of the ring, Eq. (|18|) takes the form 



In this case the upper bound is proportional to the intensity of the magnetic field, and the 
persistent current is allowed to grow indefinitely with the latter (N is now not constant). 
A persistent current of the order of ea; c /47r (we assume that the geometric parameters are 
fixed and that ln{R2/R\) is of order unity) is consistent with the results of quantitative 
calculations of the persistent current in an ideal ring [|]]. 

These examples show that the rigorous upper limit derived in this paper can provide a 
reasonably good estimate of the actual maximum value of the persistent current in a given 
toroidal system. Of course, since the sign of the current is left completely undetermined 
by our arguments, we cannot say anything about the result of averaging, for example, with 
respect to the impurity distribution in an ensemble of mesoscopic rings. 

We emphasize that the upper bound derived in this paper applies to the current in 
thermodynamic equilibrium. Our argument is not applicable to the current carried by a 
superconducting ring, because, in that case, the current- carrying state is not the true equi- 
librium state, but, rather, a metastable state which cannot decay ||. This work has been 
supported by NSF Grant No. DMR-9100988. I thank Michael Geller for useful discussions. 
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FIGURES 



FIG. 1. Cylindrical coordinates for a toroidal body. 
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